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ABSTRACT 

Observations have evidenced that passively evolving massive galaxies at high redshift are much more com- 
pact than local galaxies with the same stellar mass. We argue that the observed strong evolution in size is 
directly related to the quasar feedback, which removes huge amounts of cold gas from the central regions in a 
Salpeter time, inducing an expansion of the stellar distribution. The new equilibrium configuration, with a size 
increased by a factor > 3, is attained after ~40 dynamical times, corresponding to ~ 2 Gyr This means that 
massive galaxies observed at z > 1 will settle on the Fundamental Plane by z ^ 0.8-L In less massive galaxies 
(Mi, < 2 X 10"^ Mq), the nuclear feedback is subdominant, and the mass loss is mainly due to stellar winds. 
In this case, the mass loss timescale is longer than the dynamical time and results in adiabatic expansion that 
may increase the effective radius by a factor of up to ^ 2 in 10 Gyr, although a growth by a factor of ~ L6 
occurs within the first 0.5 Gyr. Since observations are focused on relatively old galaxies, with ages > 1 Gyr, 
the evolution for smaller galaxies is more difficult to perceive. Significant evolution of velocity dispersion is 
predicted for both small and large galaxies. 

Subject headings: Galaxies: formation - galaxies: evolution - galaxies: elliptical - galaxies: high redshift - 
quasars: general 



1. INTRODUCTION 

Several recent observational studies have found that mas- 
sive, passively evolving, galaxies at z > 1 are much more 
compact than local galaxies of analogous stellar mass (Fer- 
guson et al. 2004; Trujillo et al. 2004, 2007; Zirm et al. 2007; 
Cimatti et al. 2008; Damjanov et al. 2008). Since similarly 
superdense massive galaxies are extremely rare or absent at 
z ~ (Shen et al. 2003) a strong size evolution, by a factor 
~ 3 or more, is indicated. No convincing mechanism able 
to account for such size evolution has been proposed so far 
In the following we show that an expansion consistent with 
the observed one is naturally expected as a consequence of 
feedback from active nuclei (Silk & Rees 1998; Granato et 
al. 2001, 2004), which is now widely recognized as a crucial 
ingredient of semi-analytic models (see, e.g., Di Matteo et al. 
2008). 

In the local Universe spheroidal galaxies occupy a quite 
narrow region, the Fundamental Plane, in the 3-dimensional 
space identified by the effective or half-light radius r^, by the 
central velocity dispersion ao, and by the mean surface bright- 
ness within rg (Djorgovski & Davis 1987). The tight color- 
magnitude relation, the color-velocity dispersion relation, and 
spectral line indices imply that the bulk of stars of elliptical 
galaxies formed at z > 1.5. The enhancement of a-elements 
abundances with respect to iron in massive elliptical galaxies 
entails that most of their stars formed within the first Gyr of 
their life (see Renzini 2006 for a comprehensive discussion). 

An additional key result is the generic presence of a Super 
Massive Black Hole (SMBHs) in the center of local ellipti- 
cal galaxies. Its mass is directly proportional to the mass of 
the old stellar population, Mbh 2 x 10"^ (Magoman et 
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al 1998; see Ferrarese & Ford 2005 for a review), implying 
that quasars and spheroidal galaxies form and evolve in strict 
relation and with mutual feedback. Specifically, it has been 
suggested that the central SMBH grows until until its feed- 
back unbinds the residual gas in the host galaxy and sweeps 
it out through a high velocity wind, thus halting both the star 
formation and its own fueling and establishing a relationship 
between the SMBH mass and the stellar velocity dispersion 
(Silk & Rees 1998). Strong AGN feedback also appears to be 
the only viable mechanism to explain the exponential cut-off 
at the bright end of the galaxy luminosity function (Croton et 
al. 2006) and the observed bimodality in the color-magnitude 
diagram of galaxies at z > 1 .5 (Menci et al. 2006). Direct ob- 
servational indications of massive outflows close to high red- 
shift quasars, consistent with this scenario, have been reported 
(e.g. Simcoe et al. 2006; Prochaska & Hennawy 2008). 

As shown below, the amount of gas rapidly stripped from 
the central regions of the galactic halo can be large enough 
to drive a large increase of the galaxy size. The puffing up 
of a system by rapid mass loss is a well known phenomenon, 
extensively studied both analytically and through numerical 
simulations, with reference to galaxies (Biermann & Shapiro 
1979), and, especially, to globular clusters (Hills 1980; Good- 
win & Bastian 2006). Slower, adiabatic expansion is caused 
by mass loss due to stellar winds or supernova explosions 
(Hills 1980; Richstone & Potter 1982). 

In this Letter we first summarize the effect of mass loss 
on the galaxy size evolution (§ 2), then we present quantita- 
tive estimates on the evolution of the effective radius and of 
the central stellar velocity dispersion as a function of galac- 
tic mass, using the Granato et al. (2004) model as a reference 
(§ 3), and finally, in § 4, we summarize and discuss our results. 

2. SIZE EVOLUTION DUE TO MASS LOSS FROM VIRIALIZED 
SYSTEMS 

The effect of the mass loss on the structure and dynamics of 
a virialized stellar system depends on the amount of ejected 
mass and on the timescale of ejection. Two regimes can be 
identified, corresponding to an ejection timescale, Tej, shorter 
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or longer than the dynamical timescale Tdyn- 

2.1. Rapid mass loss 

A rapid mass loss (tej < Tdyn) results in a shallower po- 
tential well at essentially constant velocity dispersion. If M 
and M' are the initial and final masses, the final energy, E' , 
of a spherical system is related to the initial energy E by 
E' =E{M' /Mf{2-M/M') (Biermann & Shapiro 1979), so 
that if M/M' > 2 the system has positive energy and is un- 
bound. If M /M' < 2, the system will eventually relax to a 
new equilibrium configuration. If the latter is homologous to 
the initial one, the ratio of initial (R) to final (/?') radii is 

R/R' = 2-M/M'. (1) 

Numerical simulations essentially confirm this simple esti- 
mate and show that, if the system is not disrupted, the new 
equilibrium configuration is reached after about 30^0 initial 
dynamical times (see, e.g., Goodwin & Bastian 2006). 

It is important to note that the above argument assumes that 
the system is self-gravitating and isolated. Clearly the dark 
matter (DM) halo, which extends far beyond the stellar distri- 
bution, exerts a stabilizing action and prevents its disruption. 

2.2. Long lasting, adiabatic mass loss 

If the mass loss occurs on a timescale Tej > Tdyn, the system 
expands through the adiabatic invariants of the orbits of stars, 
and the expansion proceeds at a rate proportional to the mass 
loss rate (Hills, 1980; Richstone & Potter 1982). We have: 

R'/R = M/M'. (2) 

In this case no global disruption of the system is possible and 
the velocity is inversely proportional to the size. 

3. MASS LOSS AND SIZE EVOLUTION OF HIGH REDSfflFT 
GALAXIES 

Observations of (sub)-millimeter bright QSOs, thought to 
be close to the onset of the massive outflows determining 
the transition from the active star-forming phase to the un- 
obscured QSO phase, indicate very large gas masses in the 
central regions of their host galaxies (Lutz et al. 2008, and 
references therein), comparable to the stellar masses (Coppin 
et al. 2008). Similar star to gas ratios are found in the central 
regions of sub-mm bright galaxies (Tacconi et al. 2008). 

About 5% of the QSO bolometric luminosity during the 
maximum quasar activity is enough to remove the gas in about 
one Salpeter time (~ 4 x 10^ yr), through outflows of M > 
l,OOOM0yr"' (e.g. Granato et al. 2004). The high precision 
spectroscopy of the QSO SDSSJ1204 + 0221 by Prochaska 
& Hennawi (2008) is indeed consistent with a high velocity 
(y ~ l,000kms"'), massive outflow of M 3,OOOM0yr"'. 

There are thus both observational indications and theoret- 
ical arguments for the occurrence of the conditions leading 
to the expansion described in § 12.11 For sake of definiteness, 
we will present quantitative estimates exploiting the model 
by Granato et al. (2004), which reproduces the crucial obser- 
vational features of the coevolution of quasars and elliptical 
galaxies (see Lapi et al. 2006 for a comprehensive discus- 
sion). In this model, both the star formation and the growth of 
the SMBH proceed faster in more massive halos, in keeping 
with the observed "downsizing". In the most massive galax- 
ies, the quasar feedback quenches the star formation after 
~ 0.5 Gyr, thus accounting for the observed a-enhancement. 
The star formation occurs in the inner regions of the galac- 
tic halo, encompassing a fraction < 20-30% of the baryons 



associated to the halo. The star formation may be triggered 
and enhanced by the rapid merging of a few sub-halos. The 
dynamical friction is able to dissipate angular momentum of 
the collapsing baryons, transferring it to DM particles. The 
numerical simulations by El-Zant et al. (2001) showed that 
dynamical friction between baryon and DM clouds within the 
scale radius of the Navarro et al. (1997) density distribution 
(f < fs) is strong enough to overcome the expected adiabatic 
contraction of DM caused by the deepening of the potential 
well under the influence of the shrinking baryons. Eventually 
the central regions are gravitationally dominated by baryons 
(stars and cold gas), that may share the velocity field of the 
DM clouds (cfr Fig. 3 of El-Zant et al. 2001). Under more 
general assumptions, the baryon collapse is expected to result 
in stellar velocity dispersions somewhat higher than those of 
the DM. We set ct* « /o-ctdm- 

We can assume that the stars and the cold gas, of mass M^, 
andMcoid respectively, are in virial equilibrium just before the 
mass loss. Neglecting the DM contribution to the mass in 
the central region, its gravitational radius can be expressed as 
Tg « G(M* +Mcoid) /cj and the dynamical time scale of the 
collapsed baryon component is Tdyn = (7r/2)(rg/cr*). 

If we adopt a Sersic's (1968) law for the projected stel- 
lar density profile, the effective radius, that we identify with 
the projected half stellar mass radius, is related to the grav- 
itational one by = Ss{n)rg, and the observed luminosity 
weighted central velocity dispersion is related to the global 
density-weighted velocity dispersion by ctq = SKin) < cr^ >, 
where n is the Sersic index (n = 4 for a de Vaucouleurs (1959) 
r'/"* law). Both Ss(n) and Sxin) are rather weak function 
of n tabulated by Prugniel & Simien (1997); in particular, 
5,(4) ~ 0.34 and Sk(4) ~ 0.52. 

For a Navarro et al. (1997) profile, the halo circular ve- 
locity, Vh, and the virial radius rvir are well approximated 
by Vh 230(MH/10^2y^^Q)l/3[(l +z)/4]'/2kms-' and r^i, = 
8O(MH/lO'2M0)'/3[(l+z)/4]-ikpc, respectively. If the den- 
sity profile has a scale length r,, the DM velocity dispersion is 
cdm = /(c)'^^Vh, /(c) ~ 1 being a weak function of the con- 
centration c = r^u/rs- Rearranging, we obtain the dynamical 
time and the effective radius as a function of virial radius and 
mass of the host halo, and of the mass in stars and in cold gas 
at the onset of the mass loss: 

5,(4)(M,+Meold) 

The model gives the mass in stars and in gas, and the mass 
loss for any choice of the galaxy halo mass, Mh, and of the 
virialization redshift. We can then compute the effective ra- 
dius with eqs. (O and its evolution with eqs. ([T]l and (|2|l. The 
application of the virial theorem, using the coefficients 5,(4) 
and Sk{A-), allows us to derive the central line-of-sight veloc- 
ity dispersion in different evolutionary phases. The results are 
shown in Fig.|2] 

Simple analytical approximations for the mass in stars 
and cold gas as a function of halo mass and redshift were 
derived Mao et al. (2007, Appendix A) in the frame- 
work of the Granato et al. (2004) model. The key 
timescales of the model are the condensation time of the 
cold gas out of the hot phase at virial temperature, fcond ~ 
9 X lO^[(l+z)/4]-i-5(MH/lO'2M0)"-2yr, and the duration of 
the star formation phase before the onset of the quasar 
super-wind. Aw » 6 x \(fi[(\+z)m-^-''F{MH /lO^^MQ)yr, 
where Fix) = 1 for x > 1 and Fix) = jc"' for x < 1 . The mass 
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Fig. 1. — Effective radius versus stellar mass. For A/* > 2 X IO'^Mq tri- 
angles, asterisks, and squares show, respectively, the values of )>: at the end 
of the blight quasar phase; when the system has reached a new virial equi- 
librium; at the present time. For M* < 2 X IO'^Mq crosses, triangles, and 
squares show, respectively, the values of r,; just at the end of baryon collapse; 
after further 0.5 Gyr; at the present time. The thi'ee connected symbols refer 
to different halo virialization redshifts (zvk = 2, 4, and 6); smaller values of 
re correspond to higher values of Zvii- The grey area marks the itlo" local 
size-M,t relation (Shen et al. 2003). Data points as in Fig. 15 of Cimatti et al. 
(2008). A typical error bar is shown in the lower right-hand comer 



in stars formed during the burst is 

M^/Mh ~ S ^-[1 -eXp(-Afbur.st/fcond)], 



(4) 

57-1 

where ~ 0.18 is the cosmic baryon to DM density ratio and 
the function 7 ~ 1 -R + Q.6{Mh /10^^Mq)-^I\{\ +z)/4]-i in- 
corporates the supernova feedback effect. The ratio between 
the condensation time and the characteristic time of star for- 
mation s = fcond/f* ~ 5 accounts for the clumping of the gas 
distribution; R ~ 0.3 denotes the fraction of mass restituted 
by dying stars to the ISM. During the last e-folding time the 
quasar activity ejects the residual cold gas, which amounts to 



q 



[S (eXp(Afburst/fcond)- 1)] 



(5) 



For massive galaxies hosted in halos of > IO'^Mq, 
Mcoid/Mj, ^2/3. The ejection of a mass ~ Mcoid can then 
produce a factor ~ 3 increase in galaxy size (cfr eq. ([T]i). 
On the other hand, for smaller galaxies {Mh < lO'^M© and 
present day stellar mass M* < 2 x 10'*' M©) the nuclear activ- 
ity is much weaker and occurs later, when the ratio Mcoid/A^*, 
and hence Mej/M*, is small. This is in keeping with the gen- 
eral picture that the early evolution of small galaxies is ruled 
by SN feedback, while the evolution of the large ones is con- 
trolled by quasar feedback (see, e.g., Shankar et al. 2006). 

4. DISCUSSION AND CONCLUSIONS 

The variation of the effective radius implies that elliptical 
galaxies spend the initial part of their life outside the 'local' 
fundamental plane and only subsequently settle on it. 

For massive galaxies (M* > 2 x 10'"Mq), for which the 
expansion is mostly due to quasar-driven super-winds, Fig.[T] 
displays three steps in the evolution: i) the end of the quasar 
phase at cosmic time fQso — fvir + Afburst, which corresponds 
to the end of the rapid mass loss; ii) the intermediate time fint 
when the stellar component reaches a new virial equilibrium 
with a larger size; Hi) the present time to, at which a further, 
minor increase of the size is achieved due to stellar winds. 
The initial effective radius has been computed from eq. ^ 
with/,, « 1.3. 



3.0 : 

2.8 - 

2.6 - 

2.4 - 

1 2-2 : 

O 

~ 2.0 - 




log M. [Mj 

Fig. 2. — Evolution of the central line-of-sight velocity dispersion versus 
stellar mass. Symbols as in Fig.[T] the same values of Zyir and fa as in Fig.[T] 
are adopted. The grey area reproduces the distribution of local data (Gallazzi 
et al. 2006). 

Equation ^ shows that the size expansion can be extremely 
large as the fraction of expelled gas approaches 50% of the 
initial mass. However, at variance with the case of star clus- 
ters, confined by their own gravitation field, the presence of 
the DM halo, dynamically dominant at large radii, prevents 
the disruption and restricts the range of possible final struc- 
tures of large elliptical galaxies. 

The time lapse required to reach the new virial equilibrium 
can be inferred from numerical simulations of stellar clusters, 
which show that the system reaches the new equilibrium af- 
ter about 40 dynamical times (computed for the initial con- 
figuration), independently of the initial mass (Geyer & Burk- 
ert 2001; Goodwin & Bastian 2006). By scaling up this re- 
sult, massive galaxies are expected to reach a new equilib- 
rium after 2 Gyr, at fim ^ fQso + 2Gyr (asterisks in Fig.[T]i- 
Since then only minor adiabatic mass losses occur, producing 
slight changes of radius and velocity dispersion (from aster- 
isks to squares). We recall that the quasar activity statistically 
reaches its maximum at redshift z ~ 2, corresponding to a cos- 
mic time fcosm ~ 3 Gyr The rapid mass loss is then expected 
to statistically peak at the same redshift, triggering size vari- 
ations that stabilize after ^ 2 Gyr Thus the massive galaxies 
should on the average reach the local size-mass relation at 
fcosm ^ 5 Gyr, corresponding to z ~ 0.8. This result is in keep- 
ing with sizes observed at z < 1 (Cimatti et al. 2008) and with 
results of studies of the Fundamental Plane at z < 0.8 (Renzini 
2006). In the redshift interval 0.8 < z < 3 a large scatter in the 
size-stellar mass relation is expected and indeed observed. 

This scenario may appear to be contradicted by observa- 
tions of Ferguson et al. (2004), who found an average size 

~ 1.7kpc for luminous Lyman Break Galaxies (LBGs) at 
z ~ 4, while we would expect 0.5 kpc, inserting the ap- 
propriate masses in eq. (O. However as pointed out by Joung, 
Chen & Bryan (2008), the apparent effective radius could be 
a factor of about 3 larger than the intrinsic one, because of the 
presence of dust in the central star forming regions. 

In the case of rapid mass loss, the velocity dispersion of the 

-1 /2 

Stars at virial equilibrium scales as a^, (x re . As shown by 
Fig.|2] the size evolution illustrated by Fig.[T]brings the quite 
high velocities inferred from the observed values of Mj, and 
in the compact phase, to within the locally observed range. 

The structural evolution of low mass E galaxies follows a 
different path, driven by the slow mass loss due to galactic 
winds and supernova explosions, because their nuclear activ- 
ity is of low power Adopting a Chabrier (2005) Initial Mass 
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Function (IMF), we find a mass loss of a factor of ^ 2 on 
a time scale of several Gyrs (the factor is only 1 .4 for a 
Salpeter (1955) IMF); after eq. Q the size increases by the 
same factor However most of the expansion occurs when 
these galaxies are young. For a Chabrier IMF, akeady half a 
Gyr after the baryon collapse the size has increased by a fac- 
tor of 1 .6 (triangles in Fig.[T]) and the subsequent expansion 
is limited to a factor of 1.3 (squares). Thus such galaxies 
observed at ages > 0.5 Gyr should exhibit a size quite close 
to that of local spheroidal galaxies with the same stellar mass. 
As a consequence, high-z galaxies with M^, < 2 x IO^'Mq 
should also exhibit a smaller scatter in the effective radius- 
stellar mass relation. They are predicted to have significantly 
more compact sizes only at ages < 0.5 Gyr, when their star 
formation rates are of tens Mq/jt (see Fig. 1 of Lapi et al. 
2006), typical of relatively low luminosity, high-z LBGs. In 
the adiabatic expansion case the velocity dispersion scales as 
fj* (X r"'. Again, this brings the high initial velocities within 
the locally observed range. 

As apparent from Fig. [H we predict significant evolution 
for galaxies with > 2 x 10'" M0, well below the threshold 
of > 5 X 10" M0 impUed by the semi-analytic model of 
Khochfar & Silk (2006). 

In conclusion we suggest that the rapid mass loss driven by 
the quasar feedback is the main agent of the size and velocity 
dispersion evolution of massive spheroidal galaxies. Lower- 
mass galaxies experience a weaker, but non negligible evo- 
lution (of amplitude depending on the adopted IMF), due to 
mass loss mainly powered by supernova explosions; most of 
it occurs during their active star-formation phase. Although 
our calculations have been carried out in the framework of the 
Granato et al. (2004) model, this evolutionary behaviour is a 



generic property of all models featuring large mass loss due to 
quasar and/or supernova feedback. Observational evidences, 
some of which are briefly summarized in §[3] suggest that 
quasar driven high velocity outflows may have removed from 
the central regions of massive galaxies a gas mass comparable 
to the mass in stars on a timescale shorter than the dynamical 
time. If so, simple, model independent, physical arguments, 
presented in §|2] imply a swelling of the stellar distribution 
and a decrease of the stellar velocity dispersion. The model 
contributes detailed, testable, predictions on the evolution of 
the effective radius and of the velocity dispersion as a function 
of the galactic age and of the halo mass. It specifically pre- 
dicts different evolutionary histories for galaxies with present 
day stellar masses above and below = 2 x IQ^^Mq. 

The energy injected by dry mergers into the stellar systems 
can also enlarge their size. However mergers increase the 
mass in stars too, and, to first order, move galaxies roughly 
parallel to the size-mass relation (van Dokkum et al. 2008; 
Damjanov et al. 2008), while the data show size evolution at 
fixed mass in stars. 

The simple analysis presented in this Letter is intended to 
be a first exploration, sketching a promising scenario for in- 
terpreting challenging observational results. This scenario can 
be tested, on one side, by numerical simulations with appro- 
priate time resolution, properly taking into account the in- 
teractions between baryons and DM, and, on the other side, 
by observations of size and velocity dispersion especially of 
lower mass high-z galaxies, that are expected to show an evo- 
lutionary behaviour different from that of massive galaxies 
because of the different mass loss history. 

This research has been partially supported by ASI contract 
1/016/07/0 "COFIS". 
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